Towards a dynamical theory of multifractals in turbulence 
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Making use of the exact equations for structure functions, supplemented by the equations for dissi- 
pative anomaly as well as an estimate for the Lagrangian acceleration of fluid particles, we obtain 
a main result of the multifractal theory of turbulence. The central element of the theory is a dissi- 
pation cut-off that depends on the order of the structure function. An expression obtained for the 

exponents s„ in the scaling relations (ff )"/(f^)^ ^ oc Re"", between the velocity gradients §^ and 
the Reynolds number Re, agrees well with experimental data. 
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Questions of small-scale universality in fluid turbulence 
hover around the universality of the scaling exponents 
^n,o of velocity structure functions deflned through rela- 
tions such as 
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where u(x) is the velocity component along the separa- 
tion distance r, measured at the position x and e is the 
mean rate of energy dissipation. Here r lies in the inertial 
range given by 77 << r << L, where L is the large-scale 
at which the energy is being injected and rj = (z^^/e)^/^ is 
the dissipation scale, v being the fluid viscosity. The zero 
index in ^„_o shows that no powers of the transverse ve- 
locity increments are involved in this particular definition 
(1). Kolmogorov [1] assumed that the velocity fluctua- 
tions in the inertial range are independent of both L and 
T], and that e, regarded as equal to the energy flux across 
scales, is the only relevant dynamical parameter. As is 
well known, Kolmogorov's proposal yields the linear rela- 
tion ^„^o — n/Z. Since, in the limit of vanishing viscosity 
(or, as 77 ^ 0), Kolmogorov's scaling theory combines the 
exact expression [2] Sz,o — ~5^'"i if is reasonable to re- 
gard the theory loosely as dynamic. However, experimen- 
tal and numerical data in three-dimensional turbulence 
have shown (see Ref. [3] for a recent account) that the 
scaling exponents ^„_o depart from n/3, and that there 
exists a more complicated nonlinear spectrum of scaling 
exponents ^„,o- Its theoretical explanation for the veloc- 
ity field has proved to be elusive, though considerable 
progress has been made for passive scalars [4]. 

In recent past, the problem of scaling exponents in tur- 
bulence has been analyzed within a general framework 
of the theory of multifractal (MF) processes reviewed in 
Ref. [5] . This approach has led to interesting interpreta- 
tions and novel work (see [5,6] for incomplete list), but 
its shortcoming is the lack of connection with the dy- 
namical equations. In this paper, a main relation of the 
MF theory is derived from dynamical equations, supple- 
mented both by an order-of-magnitude estimate for the 



Lagrangian acceleration of a fluid particle, and the earlier 
work on dissipative anomaly [7,8] . 

For background, we review here the main ideas of the 
inertial-range MF theory, whose basis are the assump- 
tions that (a) the velocity increments 5rU have the form 



(5ru(x) ?i(x + r)— u(x) r, 
— oc (- 



(2) 



where u' ^ 6lu may be regarded as the root-mean-square 
value of u, and (b) there exists a spectrum of exponents 
h related to the fractal dimension of their support D{h). 
Thus, [r/LY^^'^'^'^ is proportional to the probability of 
the velocity increment falling within a sphere of radius r 
on a set of dimension D{h). It is clear from (2) that 
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where dfi(h) is the weight of a local value of exponent 
h. Thus the scaling exponents ^n,o are directly related to 
the spectrum D{h). The goal of the theory is to find the 
functions D{h) and fJ,{h). If one evaluates the integral in 
(3) in the steepest descent approximation, as in the stan- 
dard procedure, the precise form of /i(/i) is irrelevant and 
only the spectrum D{h) needs to be determined. How- 
ever, this cannot be done within the MF theory itself. 

Multifractality in the inertial range will have conse- 
quences for dissipation scales as well. The authors of 
Ref. [9] used the local scaling (2) to construct eddy- 
turnover times that depended on h, equated them to dif- 
fusion times scales rf /v, and showed that a spectrum of 
/i-dependent dissipation scales can be written as 
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where the large-scale Reynolds number Re ~ u' L/v. The 
exponents in (4) have to be related somehow to the spec- 
trum of scaling exponents of structure functions. This 
can be done by assuming, for any h, that (2) is valid 
only for scales r > rj{h), with smoothness for smaller 
scales [10,5]. One can then evaluate (3) for scales larger 



than r]{h), using the steepest descent approximation up 
to the cut-ofF scale r = ri{h). It is easy to show [5] that 
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where Cn,o = pin, 0) — ^ and p{n, 0) is the solution of 
p{n,0) = 2n — £,p,o. Our specific goal is to obtain s„ 
theoretically. 

A brief remark on our strategy may be helpful here. 
If the structure functions Sn,m = {Sru)"{drv)™, where v 
now is the velocity component normal to the displace- 
ment vector r, of the form An^mT^"''"^ , are non-analytic 
in the inertial range, and the viscous dissipation is the 
only mechanism for smoothing the singular nature of the 
structure functions in the inertial range, the balance be- 
tween them occurs at the length scale r — > ?7„.o, where 
rin_o = rjn is an order-dependent length scale nominally 
separating the analytic and singular intervals. The ana- 
lyticity of structure functions in the viscous range yields 
Sn.o Of (9a;u(0))"r", so we have [11] 
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This equation defines the field 77(x, t) through moments 
of velocity gradients, and picks out the strongest singu- 
larity of a chosen order dominating the inertial range 
asymptotics. Our strategy is based on the idea that if 
an n-th order structure function evaluated at the appro- 
priate cut-off is Sn,oiiln,o) = ^n??^"'", with the Reynolds- 
number- independent proportionality coefhcient An , then 
Sn.oi^) '^ r^"'" for r in the inertial range. (Henceforth, to 
simplify notation, we will often set e = L = 1 and omit 
the second index in the ^'s and ??'s.) 

As the first step in the theory, we write the exact equa- 
tion for structure function of order 2n [12,11] (see also 
Refs. [13,3,14]) as 

dS2n.o d-l^ _ (d-l)(2n-l)^ 

7[ \ 'J2ri,0 — J2n-2,2 

or r r 



+ {2n - l)(5^a^C/2"-2, (7) 

where the increment of the r-component of Lagrangian 
acceleration of a fluid particle is given by 

Sra^ - -[dMx')-dMx)] + ''[^lM^')-^K^)] (8) 

and x' = X + r. 

The second step requires the closure of Eq. (7), for 
which we need an expression for SrCix in terms of veloc- 
ity increments. The Laplacian in (8) can be represented 
in terms of finite differences on the dissipation cut off 77 
and, by virtue of Eq. (6), the acceleration increment can 
be made a function of two fluctuating variables (opera- 
tors) 6riU and rj. To make further progress, however, it is 
necessary to express 77 in terms of the velocity field itself. 

We consider two scenarios. In the first, we have the 
option of expressing the acceleration terms through ei- 
ther a model for the conditional mean of the acceleration 



increment for a fixed value of the velocity increment, or 
through a direct relation between Sa and 6u — somewhat 
in the spirit of Kolmogorov's refined similarity hypoth- 
esis [15]. Choosing the former option, we model in the 
limit r —> rj where 77 is a generic local dissipation scale, 
the x-component of acceleration term as [16] 
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where r « rj/SriU is the life-time of a fluctuation on the 
scale 77. In the last step in Eq. (9), we have used 
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(10) 



where rj is to be regarded as a random field. This step 
reduces the number of random fields from 2 to 1. Ex- 
pression (10), which is central for our theory, is proposed 
here on dimensional grounds but will be obtained below 
from a second scenario considering dissipative anomaly. 
This second scenario follows Polyakov's work [7] on sta- 
tistically steady turbulence due to the one-dimensional 
Burgers equation stirred by a large-scale random force. 
In that work, on the basis of the energy balance equation 
for 7^ = 0, namely. 
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Polyakov derived the dissipation anomaly, which is re- 
lated to the local form of the Kolmogorov law [2] as 
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0, and D^ -F = /(x+f )u(x-f )+/(a;- 

0, while D — v(u(x + y)d^_yu{x — 

y)dl^yu{x + y)) when the forcing / is zero. 



Here, y ^ rj 
^)u{x + |) when ly 
y) + u(x - 

Equation (12) balances the singular contributions in the 
limit 7/ ^ 77 ^ 0, while the regular contributions disap- 
pear by virtue of (11). The coordinate shift y in Eq. (12) 
is identical to Kolmogorov's displacement r = y ^ rj. 
If, as 7; ^ 0, the velocity field is non-differentiable (i.e., 
singular), the left side of Eq. (12) does not approach zero 
even in the limit v ^ 0. In a statistically steady state, 
Eq. (12) immediately gives (u{x + y) — u{x)Y — —12Fy 
for the inertial range y ~ r << L. We can see that the 
celebrated — 4/5-ths law of Kolmogorov [2] is not locally 
valid because of the 0{dxU^) term in (12); this term can, 
however, be eliminated by averaging (12) over the direc- 
tions of velocity vector u/u. 

Now using the finite difference definition of all deriva- 
tives on a dissipation scale 77 — > 0, we can write, after 
some algebra, that 

1 ■u^(x+) - 7/^(x") (?i(a:; + 277) - 7i(x))^ 
3 77 677 

« [2i^/f]'^] X [7i^(x+) + u'^{x^) + u{x+)u{x - 377) 

+u{x^)u{x + 3r])-4:u{x+)u{x^)], (13) 



where x^ = x zt rj. This equation is correct up to 
O(rj^). As mentioned ear her, the single-point contribu- 
tion to this relation disappears when averaged over the 
"directions" of rj. While the left side of (13) involves 
two-point differences, the right side includes contribu- 
tions from four shifted points. To proceed further, we 
assume that 77 plays the same role as the width of typical 
shock structures, and conclude that u{x + 877) — u{x) « 
u{x^) — u{x^) « u(x + 2r]) — u{x) , as a result of which the 
right side of (13) is 0[t/((5^u)V??^]- This leads to (10). 

In three dimensions, however, additional terms appear 
due to the pressure gradient- velocity product. The rele- 
vant extensions have been made in Refs. [8]. The finite- 
difference representation of the equations from Refs. [8] 
on the dissipation scale r/ — > yield the estimate 
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(14) 



Since on the dissipation scale the pressure and dissipa- 
tion contributions are of the same order [16], expression 
(14) gives the same balance relation, v « r]SriU, obtained 
above. Thus, the relation (10) applies also to three- 
dimensional turbulence. 

Two comments are in order. First, the model for ac- 
celeration should include the 0{{Sriu)'^ /rj) quadratic con- 
tribution coming from the pressure terms [16]. However, 
the pressure term simply renormalizes the coefficients in 
front of the remaining contributions to (7) and, as a re- 
sult, does not alter the steps presented above (see also 
Ref. [11]). Second, to our knowledge, there are no ex- 
perimental or numerical data that directly address the 
conditional acceleration term of Eq. (9), though related 
conditional data are accumulating rapidly [17,18]. 

Substituting Eq. (9) (after using Eq. (10)) into Eqs. 
(7) and (8) we obtain an infinite set of equations cou- 
pling the structure functions 52n('') and S2n+i{r). These 
equations are valid for all magnitudes of displacement 
r ^ L, including r -^ rj2n- It will become clear below 
that 772ri ^ iT2n+i and, as a result both S2n{'>]2n) and 
'S'2ri+i('72n) are in their respective algebraic ranges, i.e., 
S2n{mn) oc ?7^^" and S'2„+i(?72n) oc ??2n""^'- Thus, on the 
scale ?72n, Eqs. (9) and (7) give 



4l 



We thus have 
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(15) 



(16) 



giving, for n = 1, the well-known relation [19] for the 

1 
dissipation scale 772 ~ i?e«2-2. Since by Holder inequal- 
ity, for all q > p, the exponents £,p > £,q, it follows from 
(16) that ?72n > rj2n+i, which justifies the derivation of 
expression (15). By virtue of Eq. (6) we have 
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(17) 



exponent 


S2 


S3 


S4 


S5 


S6 


theory 





0.066 


0.176 


0.34 


0.49 


experiment 






0.05 
(0.06) 


0.20 
(0.16) 


0.32 
(0.31) 


0.54 
(0.50) 



TABLE I: The scaling exponents s„ from theory and experi- 
ment. The numbers in parentheses are from Ref. [18]. 



To compare Eq. (17) with the outcome of the multifractal 
formula (5), we notice that both are the same in the 
limit n ^ 1. In the limit n ^ 00, if the exponents 
S,n -^ Coo = const, or C„ -^ an, Eq. (17) and the outcome 
of mutifractal formula (5) are identical. 

The moments of velocity derivatives are given from Eq. 
(17) to be 



^n = (|^)"/(|^)^'oci?e-. 
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S2r. 



6 



2n 



S,2n — S,2n+1 — 1 



(18) 



(19) 



This expression can be evaluated readily if we know ^„ 
Using the result obtained in [11,12] for C„, we get 



Cn 



1.15n 



3(H-0.05n)' 



(20) 



The exponents S2-sq from (19), after using (20), are 
listed in Table I. Almost the same results are obtained if 
experimental values [3] are chosen for C„ instead of (20). 
These values also agree very well with results from several 
phenomenological MF models. 

Several experimental measurements of s„ are available 
in the literature. For an incomplete list, see Ref. [18, 
20, 21]. We compare in Table I the theoretical numbers 
above with the data of [21] in the atmospheric boundary 
layer at very high Reynolds numbers and the latest wind 
tunnel measurements [18] in grid turbulence. The differ- 
ences between the two sets of experimental numbers are 
a measure of uncertainty in the data. Keeping this in 
mind, we may regard the agreement with the theoretical 
values to be very good. The conclusions from other data 
sets [20] are quite similar. 

As an aside, we note that the theory can also be used to 
show that the results of Kolmogorov's Refined Similarity 
Hypotheses (RSH) [15] are at least numerically close to 
the ones derived above. We can evaluate the moments 
of velocity derivative {dxu)"^"- by extrapolating RSH to 
rj. The dissipation averaged on this scale is of the same 
order as the unaveraged dissipation. This assumption, 
commonly adopted in the literature, appears reasonable 
if there is no structure for scales smaller than rj. From 



dissipation anomaly, we obtain 



e oc 



where 
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(22) 



is the dissipation rate averaged over a "ball" of a radius r] 
with the center at x. Equation (21) is an order of magni- 
tude estimate averaged over the "universal" Kolmogorov 
noise V = ?7e^/((5,,u)^. Combining Eqs. (21) and (18) 
for the dissipation scale of the 3n-th moment of velocity 
difference we obtain 



instead of Eq. (19), where 
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As we see, the two sets of formulae (18)-(19) and (23)- 
(24) are identical when n = 1 leading to the exact rela- 
tion {dxu)^ oc Re. With the relation (20) for the scal- 
ing exponents, both relations have the same asymptotics 
G2n — *■ i?e^" in the limit n —> oo. In the interval n > 1, 
the formulae (18)-(19) and (23)-(24) differ by no more 
than a few percent. 

In summary, the theory developed here combines the 
exact equations (7) and (8) with relations (9) and (10). 
Together, they lead to Eqs. (18)-(19), which form the 
main result of the paper. While this form is known from 
the MF theory, the present paper obtains the exponents 
theoretically and the results agree well with experiments. 

It is useful to restate here the approximations involved 
in derivation of (9) and (10). Expression (10) was ob- 
tained through the model (9) for the acceleration terms 
in (7), and also as an order-of-magnitude estimate from 
the equations for dissipation anomaly. Since, at the dis- 
sipation scale Tj, the pressure contribution simply renor- 
malizes the coefficients in the left side of equation (7), 
the expression (9) for the viscous friction force introduces 
i?e-dependence into (7). This is the reason why, for the 
fixed magnitude of the inertial range displacement r, the 
structure functions Sn,m(j) are i?e-independent, while 
the moments of derivatives are strongly i?e-dependent. 
Though we cannot prove that either scenario leading to 
(10) is rigorous, the good agreement observed with ex- 
perimental data gives us some confidence that the theory 
is a step in the right direction. 

We thank A. Bershadskii, T. Gotoh and A.M. Polyakov 
for helpful discussions. 
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